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摘 要 : 土地 利用 空间 格局 变化 反映 土地 利用 合理 性 和 生态 环境 情况 。 以 位 于 农 牧 业 交 错 带 的 生 


态 环 境 脆弱 区 陕 北 榆林 市 为 研究 区 ,运用 核 密度 分 析 法 表达 土地 利用 空间 格局 变化 ,采用 土地 利 
用 动态 度 、 土 地 利用 综合 动态 度 、 土 地 利用 变化 强度 、 土 地 利用 转移 矩阵 指数 分 析 榆 林 市 近 40a 的 
土地 利用 变化 速度 、 强 度 、 结 构 及 其 影响 因素 。 结 果 表 明 : 榆 林 市 近 40 a 土地 利用 变化 以 耕地 面积 
减少 建设 用 地 大 幅度 增加 和 和 未 利用 地 减少 为 主 。 东 南部 耕地 密度 等 级 降低 ,西南 部 呈现 相反 的 
变化 趋势 ;林地 呈 “ 多 核心 增长 ;草地 高 密度 核心 位 于 北部 ,等 级 升 高 ;建设 用 地 高 密度 核心 位 于 


西 


9 北部 6 个 县 区 ,增长 幅度 大 。 未 利用 地 高 密度 核心 为 西北 部 毛乌素 沙漠 , 斑 块 破碎 化 , 沙 地 治理 


效果 显著 。 政 策 因 素 、 社 会 经 济 因素 和 人 口 因素 的 综合 作用 是 影响 土地 利用 空间 格局 变化 的 主要 
原因 。 退 耕 还 林 还 草 使 得 人 地 矛盾 得 到 改善 ,减缓 了 水 土 流失 ;毛乌素 沙 地 治理 成 效 显著 。 未 来 
发 展 需 注 重生 态 格 局 安全 ,以 生态 保护 为 前 提 , 实 现 科学 绿色 、 可 持续 的 高 质量 发 展 。 
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土地 是 人 类 生活 的 物质 基础 ,土地 利用 变化 是 
人 类 活动 对 地 表 自 然 生 态 系统 影响 的 最 直接 信和 号， 
是 人 类 社会 经 济 活动 与 自然 环境 连接 的 纽带 ,人 与 
自然 的 矛盾 是 生态 和 环境 问题 的 根源 。 工 业 化 、 
城市 化 的 增长 与 资源 配置 的 增强 对 区 域 土地 利用 
空间 结构 产生 了 深远 影响 ,不 合理 的 土地 利用 
变化 造成 景观 破碎 化 ,对 生态 系统 退化 产生 负面 影 
响 ”。 土 地 利用 变化 时 空 特征 、 变 化 的 驱动 力 机 制 、 
变化 模型 模拟 以 及 生态 环境 效益 问题 已 经 成 为 研 
究 焦点 “”。 研 究 不 同时 期 土地 利用 变化 特征 , 揭 
示 土 地 利用 变化 内 在 规律 和 影响 因素 ,以 便 合 理 优 
化 配置 土地 利用 结构 。 

榆林 市 处 于 北方 农 牧 交错 带 , 生 产 方式 和 植被 
类 型 等 具有 过 渡 性 特点 ,属于 典型 生态 脆弱 区 。 农 
牧 交错 带 土 地 利用 类 型 复杂 ,是 生态 安全 的 重要 屏 
障 , 地 理 位 置 和 界线 不 断 变 迁 ,变迁 区 内 土地 利用 
和 景观 格局 变化 复杂 ”。 由 于 能 源 经 济 产业 化 和 
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生态 恢复 过 程 双 重 作 用 ,土地 生态 安全 问题 越 来 越 
突出 ”。 土 地 利用 空间 格局 变化 直接 影响 生态 环 
境 稳定 和 安全 ”"。 前 人 基于 榆林 市 土地 利用 的 研 
究 涉 及 归 一 化 植被 指数 (NDVI) 时 空 变化 ”土地 生 
态 风 险 ”土地 承载 力 ” .生态 环境 效应 ”土地 利 
用 变化 驱动 力 ” 等 ,得 出 榆林 市 生态 环境 整体 呈 改 
善 趋势 ,但 存在 显著 区 域 差异 ,人 研究 结果 集中 于 土 
地 利用 产生 的 生态 响应 结果 ,目前 缺少 土地 利用 空 
间 格 局 分 布 差异 性 以 及 长 时 间 尺 度 土地 利用 演变 
趋势 的 研究 。 各 土地 利用 类 型 变化 产生 的 生态 效 
益 各 不 相同 ,因此 需要 基于 土地 利用 空间 格局 探究 
生态 效益 变化 背后 土地 利用 的 转变 以 及 影响 因 
素 。 核 密度 分 析 法 利用 数据 的 空间 属性 研究 空间 
数据 分 布 特征 ,是 空间 分 析 中 运用 最 为 广泛 的 非 参 
数 估 计 技 术 ; 将 地 类 属性 以 点 要 素 形 式 计 算 ,根据 
距离 越 近 权 重 越 大 的 原则 将 点 要 素 生成 平滑 的 等 
值 线 ,直观 的 反映 地 类 空间 扩散 规律 ; 核 密度 值 越 
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榆林 市 自 1956 年 开始 治理 沙 地 ,兴建 防风 林 
带 , 引 水 拉 沙 , 引 洪 淤 地 ,毛乌素 沙 地 实施 大 规模 土 
地 整治 。 退 耕 还 林 还 草 工 程 在 1999 年 发 起 以 期 减 
少 水 土 流失 并 改善 黄土 高 原 地 区 的 生态 条 件 ,是 
我 国 最 大 的 土地 退耕 造林 计划 ,荒山 造林 , 封 山 育 
林 , 封 山 禁 牧 ,促进 生态 自然 恢复 。 因 此 ,本 研究 结 
合 榆林 市 地 域 特殊 性 ,利用 核 密 度 分 析 法 ,分 析 近 


陆 性 季风 气候 ,四 季 分 明 ,日 差 较 大 ,无霜期 短 ,年 
平均 气温 8 % ,年 平均 降水 400 mm, 日照 充足 ERE 
资源 丰富 ,年 日 照 时 数 2600~2900 h, 具 有 春季 干燥 
多 风 夏季 炎热 短促 、 秋 季 多 暴雨 .冬季 干 冷 漫长 的 
特点 。 研 究 区 地 人 处 毛乌素 沙 地 和 黄土 高 原 的 过 渡 
带 , 内 蒙古 高 原 和 黄土 高 原 的 过 渡 区 ,生态 环境 脆 
弱 。 传 统 产 业 以 农 牧 业 为 主 ,能 源 矿产 资源 富 集 ， 
已 发 现 8 大 类 48 种 矿产 ,以 煤 、 气 、 油 、 盐 最 为 丰富 ， 


40 a 的 土地 利用 空间 格局 变化 ,自然 和 社会 影响 因 
素 ,揭示 土地 利用 变化 与 影响 因素 间 的 关系 ,以 期 
为 我 国 北方 农 牧 交 错 傍 生态 脆弱 区 高 质量 发 展 提 
供 科 学 依据 。 


1 研究 区 概况 及 数据 来 源 


1.1 研究 区 概况 

榆林 市 位 于 陕西 省 最 北部 ,地 理 范围 在 36°49'~ 
39°35' N, 107°15'~111°15'E Z [al ,面积 约 42921.1 
km’, 人口 为 383.84x10' 人 (图 1)。 市 辖 1 市 2 区 9 
县 ,分 别 为 神木 市 . 榆 阳 区 .横山 区 . 府 谷 县 . 靖 边 
H ELE REA KRE AER REE 清润 县 、 
子 洲 县 。 地 貌 大 体 以 长 城 为 界 ,北部 为 风沙 草滩 
区 , 占 总 面积 的 42% ,南部 为 黄土 丘陵 沟 冤 区 , 占 总 
面积 的 58% ,地 势 西高 东 低 ,境内 平均 高 程 为 
1218.63 m。 梁 状 低 山 丘 陵 区 主要 分 布 在 西南 部 山 
区 一 带 无 定 河 ` 大理 河 、 延 河 治 河 的 发 源 地 , 占 全 
市 面积 11.55%。 气 候 属 暧 温带 至 中 温带 半 干 旱 大 
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图 1 榆林 市 行政 区 划 及 地 势 分 布 图 
Fig. 1 Administrative division and geographical 
distribution of Yulin City 


目前 已 经 发 展 成 为 我 国 重 要 的 能 源 化 工 基地 。 榆 
林 市 水 土 流失 严重 ,部 分 区 域 环 境 质 量 差 ,改善 难 
度 大 ,生态 环境 脆弱 ,发 展 质量 有 待 提高 。 
1.2 数据 及 来 源 

选 1980 年 (家 庭 联 产 承包 责任 制 .沙漠 治理 成 
效 显 现 ).2000 年 (退耕 还 林 还 草 实施 ).2018 年 ( 退 
耕 还 林 还 草 基本 完成 .生态 保护 新 起 点 ) 为 时 间 节 
点 ,分 析 榆 林 市 近 40 a 土地 利用 时 空格 局 变化 。 土 
地 利用 数据 来 源 于 中 国 科学 院 资 源 环 境 科学 数据 
中 心 (http://www.resdc.cn/)。NDVI 数 据 来 源 于 (http: 
/Iwww.resde.cn) ,三 次 产业 数据 和 人 口 数据 来 源 于 
陕西 省 统计 局 官网 (http://tjj.shaanxi.gov.cn/)。 植 被 
覆盖 度 计 算 公 式 为 :f=(NDVI- NDVLYNDVL。- 
NDVInin) 0 


2 研究 方法 


21 核 密 度 分 析 法 

核 密度 分 析 法 是 计算 要 素 在 其 周边 邻 域 中 密 
度 的 非 参 数 舍 计 空间 分 析 方法 ””。 计 算 每 个 输出 
栅 格 像 元 周围 点 要 素 密 度 , 由 边界 参数 、 点 权重 和 
核 函数 确定 小 曲面 ,每 个 点 对 应 小 曲面 ,重生 小 曲 
面 数值 累加 形成 大 曲面 ,大 曲面 数值 归 一 化 生成 核 
密度 ,在 空间 上 显示 数据 点 的 聚 类 情况 。 采 用 Arc- 
GIS 10.6 空间 分 析 工 具 核 密度 ,计算 具有 土地 利用 
数据 属性 的 点 密度 。 划 分 地 类 区 域内 相对 密度 等 
级 ,反映 用 地 类 型 在 空间 尺度 的 集中 程度 和 时 间 尺 
度 的 变化 趋势 汪 。 核 密度 等 级 越 高 ,土地 利用 空间 
分 布 密度 越 大 。 


i ee 
Ae) ASKS") a) 


式 中 :ftw) 为 空间 位 置 x 人 处 的 核 密度 计算 函数 ;n 为 分 
析 范 围 内 的 点 数 沁 为 范围 国 值 ;K 为 默认 的 权重 核 
函数 ;x-x; 为 点 x 到 点 x 之 间 的 距离 。 方 程 空间 的 几 
何 含义 为 密度 值 在 每 个 分 析 窗 口中 心 点 处 向 外 递 
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减 , 当 距 离 中 心 达到 一 定 阐 值 范围 函数 值 为 0”。 
2.2 土地 利用 变化 指数 

(1) 单一 土地 利用 动态 度 指 在 一 定时 间 段 内 茶 
种 土地 利用 类 型 的 数量 变化 幅度 和 速度 ”。 地 类 
面积 变化 率 表达 增 减 幅度 和 速度 ,不同 时段 动态 度 
对 比 反 映 变 化 剧烈 程度 。 单 一 土地 利用 动态 度 指 
数 (K) 计 算 公 式 如 下 : 

U,-U, i 
k=! 7 xx 100% (2) 


式 中 :1 到 分 别 为 研究 时 段 内 期 初 和 期 末 某 一 用 地 
类 型 的 面积 ;7 为 研究 时 段 时 长 (a)。 
综合 土地 利用 动态 度 指数 计算 公式 如下: 
LC= YALU YLU,xEx100% (3) 


式 中 :ALU, 为 研究 期 内 第 a 种 用 地 类 型 转变 为 非 a 
类 用 地 类 型 的 绝对 数量 ;LU, 为 研究 期 初 第 a 类 用 地 
类 型 的 数量 ;n 为 地 类 数量 ,n=6;7T 为 研究 期 时 长 
(a). 

(2) 土地 利用 变化 强度 指数 测度 不 同类 型 土地 
变化 的 强 弱 程度 ”中 。 指 某 一 空间 单元 内 ,土地 利用 
类 型 在 单位 面积 上 发 生 的 变化 ,是 对 该 空间 单元 土 
地 面积 变化 速度 的 标准 化 ”。 

S=U,x100/UxT7) (4) 
式 中 :5 为 某 一 类 型 土地 利用 变化 的 强度 ,其 值 越 大 
表明 该 种 用 地 形态 变化 越剧 烈 ;0; 为 第 i 种 土地 利用 
类 型 在 研究 时 段 内 的 绝对 变化 量 ;U 为 各 类 用 地 类 
型 的 变化 总 量 ;7 为 研究 期 时 长 (a)。 

(3) 土地 利用 转移 矩阵 反映 某 土地 类 型 在 初期 
和 末期 与 各 土地 类 型 面积 之 间 相 互 转变 的 动态 过 
程 信息 ,包含 区 域 某 时 间 点 的 各 土地 类 型 面积 信 
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息 ,还 有 各 土地 类 型 在 期 初 和 期 未 转 出 和 转 入 的 面积 
信息 ”。 基 于 ENVI 变 化 检测 得 出 的 各 时 期 用 地 类 型 
转移 矩阵 ,制作 土地 利用 转移 矩阵 ,分 析 地 类 间 的 面 
积 转移 单一 地 类 转 出 转 和 人 和 地 类 净 变 化 等 信息 。 


3 ”结果 与 分 析 


3.1 土地 利用 空间 分 布 

表 1 和 图 2 表明 ,耕地 分 布 广 ,2018 年 耕地 面积 
占 总 面积 36.95% , 较 1980 年 降低 2.09%。 林 地 面积 
占 比 较 少 ,空间 分 布 集中 于 东南 部 ,其 他 区 域 零散 
分 布 ;2018 年 林地 面积 比例 为 5.53% , 较 1980 年 增 
长 1.24%。 草 地 密度 分 布 最 广 , 以 东北 部 最 显著 , 面 
积 比 例 最 大 ,2018 年 草地 面积 为 1885373.64 hm’, 占 
43.91% , BE 1980 年 增加 2.56%。 水 域 面积 比例 较 
小 。 建 设 用 地 集中 于 各 县 区 中 心地 带 ,1980 .2000、 
2018 年 建设 用 地 面积 比例 分 别 为 0.33% 、0.38%、 
2.55% 。 未 利用 地 主要 是 西北 部 毛乌素 沙 地 ,面积 
居 第 三 位 ,2018 面积 较 1980 年 减少 3.85%。 
3.2 土地 利用 空间 格局 变化 分 析 

图 3a 表 明 ,耕地 核心 高 密度 集中 于 西南 部 靖 边 


表 1 1980—2018 年 土地 利用 面积 
Tab.1 Area ofland use during 1980 一 2018 /hm 


土地 利用 类 型 1980 年 2000 年 2018 年 
耕地 1676427.78 1670469.43 1586294.27 
林地 184029.21 200908.62 237310.11 
草地 1776112.64 1888808.55 1885373.64 
水 域 54891.90 53572.59 51292.62 
建设 用 地 14314.32 16358.67 109622.52 
未 利用 地 589304.32 464692.33 423749.17 


À (c) 20184 
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图 2 1980 一 2018 年 土地 利用 空间 分 布 
Fig.2 Spatial distribution of land use during 1980 一 2018 
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(a) 耕地 核 密度 分 布 
1980 年 


(d) 建设 用 地 核 密度 分 布 
1980 年 N 
人 
0 70km 
(e) 未 利用 地 核 密度 分 布 
1980 年 N 
人 


0 70 km 0 70 km 
1 低 密度 区 m 较 低 密度 区 mm 中 密度 区 am 较 高 密度 区 am 高 密度 区 


图 3 土地 利用 核 密度 分 布 
Fig.3 Kernel density distribution of land use 
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县 ` 定 边 县 和 东南 部 大 部 分 县 区 ,西北 部 沙漠 地 区 
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速度 最 快 ,但 是 2000 年 后 减少 幅度 处 于 较 高 水 平 ， 


耕地 面积 较 少 ;东南 部 耕地 高 密度 与 较 高 密度 


居 第 二 。1980 一 2000 年 和 2000 一 2018 年 2 个 阶段 


2000 一 2018 年 呈 缩 减 趋势 ,高 值 密度 向 低 值 密度 等 
级 转化 ;西南 部 沙漠 区 呈现 相反 的 变化 趋势 , 较 高 
密度 和 高 密度 等 级 耕地 面积 均 增长 。 图 3b 表 明 , 林 
地 高 密度 主要 集中 于 东南 部 ;1980 一 2000 年 东南 部 
林地 面积 小 幅度 增加 ;2000 一 2018 年 ,林地 密度 核 
心 增多 ,东南 部 密度 等 级 提高 , 较 高 密度 与 高 密度 
密度 面积 增 大 。 图 3c 表 明 , 草 地 密度 核心 位 于 北部 
和 西南 部 ,以 北部 核心 为 主 ; 1980 一 2000 年 ,北部 核 
心 西 南 扩展 ,西北 部 较 低 密度 向 中 密度 等 级 提升 ; 
2000—2018 年 ,西北 部 草地 面积 增 大 , 沙 地 治理 成 
效 显 著 , 东 北部 草地 密度 增 大 , 较 高 密度 等 级 向 高 
密度 等 级 转移 ;西南 部 区 域 草 地 等 级 逆向 变化 与 耕 
地 面积 增长 有 关 。 图 3d 表明,1980 年 建设 用 地 密度 
等 级 低 ;2000 年 榆 阳 区 建设 用 地 上 升 为 高 密度 ; 
2018 年 西北 部 6 个 县 区 密度 大 幅度 升 高 ,高 密度 面 
积 大 范围 增加 ; 榆 阳 区 .神木 市 . 府 谷 县 为 高 密度 分 


综合 动态 度 分 别 为 0.14% 和 0.44% ,表明 1980 一 
2000 年 土地 利用 变化 较 平 稳 , 土 地 利用 格局 较 稳 
定 ,2000 一 2018 年 变化 较 剧 烈 ,综合 利用 程度 升 高 。 

3.3.2 土地 利用 变化 强度 ”1980 一 2000 年 草地 变化 
强度 最 大 ,其 次 是 耕地 林地、 未 利用 地 。2000 一 
2018 年 仍 是 草地 变化 强度 最 大 ,其 次 是 耕地 ET 
用 地 、 林 地。2000 年 前 ,耕地 、 林 地、 草地 为 土地 利 
用 变化 结构 的 主体 ,2000 年 后 ,耕地 变化 强度 升 高 
1.05% ,林地 变化 强度 升 高 ,草地 变化 强度 减弱 ,建设 
用 地 变化 强度 升 高 0.74% ,未 利用 地 变化 强度 减弱 ， 
耕地 .草地 建设 用 地 成 为 土地 利用 变化 结构 的 主体 。 
3.3.3 土地 利用 转移 矩阵 ”1980 一 2000 年 ,耕地 向 
草地 转移 1.63% ;林地 向 草地 .耕地 和 未 利用 地 分 别 
转移 8.22% 、1.87% 、1.31% ; 草地 向 耕地 转移 
30760.72 hm’, 占 1.63% ;未 利用 地 向 草地 转移 
2.08% , 占 比 最 大 ( 表 3)。 上 述 转移 过 程 显示 :耕地 


布 核心 区 。 图 3e 表 明 ,1980 年 未 利用 地 高 密度 分 布 
成 片 联结 ;2000 一 2018 年 沙 地 较 高 等 级 面积 大 幅度 
缩减 ,高 密度 斑 块 分 散 , 沙 地 治理 从 边缘 向 核心 地 
带 渗透 ,效果 显著 。 

3.3 土地 利用 变化 指数 分 析 

3.3.1 土地 利用 动态 度 “1980 一 2018 年 ,耕地 水域 
和 未 利用 地 减少 , 呈 负 动态 度 ; 林 地 、 建 设 用 地 增 
加 ,呈正 动态 度 ;草地 先 增 后 减 ,动态 度 呈 现 先 正 后 
负 ( 表 2)。 耕 地 2000 一 2018 年 动态 度 大 ,变化 剧 
烈 。 林 地 2000—2018 年 动态 度 高 ,变化 剧烈 ,增长 
幅度 大 ,是 1980 一 2000 年 的 2 倍 。2000 年 后 草地 负 
动态 度 0.01% ,减少 幅度 较 小 。2000 一 2018 年 建设 
用 地 达到 31.67% , 居 第 一 ,显示 该 阶段 建设 用 地 增 
长 幅度 大 。 未 利用 地 持续 减少 ,1980 一 2000 年 减少 


需求 促使 部 分 林 草 地 向 耕地 转移 ;林地 破坏 比较 严 
重 ; 沙漠 治理 使 得 未 利用 地 向 草地 恢复 。2000 一 
2018 年 ,耕地 向 草地 和 未 利用 地 转 10.96% 1.33%, 
向 草地 转移 面积 比例 最 大 ;林地 向 耕地 和 草地 分 别 
转移 17.08%F 9.54% ;草地 向 耕地 转移 11.88% ;未 
利用 地 向 草地 和 耕地 分 别 转移 5.15% 和 1.09%。 上 
述 转移 过 程 显示 :退耕 促使 耕地 向 林 草 地 转移 , 且 
效果 显著 ;存在 林 草 地 破坏 现象 ;未 利用 地 恢复 成 
效 明 显 。 

3.4 土地 利用 变化 影响 因素 分 析 

3.4.1 和 气候 因素 ”1980 一 2018 年 平均 气温 呈 增 长 趋 
势 , 降 水 量 波动 较 大 。 榆 林 市 地 处 干旱 区 ,对 气温 
和 降水 变化 较 敏 感 。 农 牧 交错 带 以 降水 量 .干燥 度 
等 为 依据 3”, 其 范围 正在 发 生 时 空 变化 *, 显 示 气 


表 2 土地 利用 变化 指数 
Tab.2 Index of land use change 


1980 一 2000 年 2000 一 2018 年 
土地 利用 类 型 - — - 
面积 变化 /hm 单一 动态 度 /% 变化 强度 /% 押 积 变化 /hm 单一 动态 度 /% 变化 强度 /% 
耕地 -5958.35 -0.02 0.70 -84175.20 -0.28 1.75 
林地 16879.41 0.46 0.47 36401.48 1.01 0.53 
草地 112695.90 0.32 3.44 -3434.91 -0.01 2.17 
水 域 -1319.31 -0.12 0.03 -2279.97 —0.24 0.07 
建设 用 地 2044.35 0.71 0.05 93263.84 31.67 0.79 
未 利用 地 -124612.00 -1.06 0.31 -40943.20 -0.49 0.24 


202108.00028v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


tió wn 


表 3 土地 利用 转移 矩阵 


Tab.3 Land use transfer matrix 


土地 利 1980 一 2000 年 2000 一 2018 年 

用 类 型 。 耕地 林地 草地 水 域 ”建设 用 地 未 利用 地 ”耕地 林地 草地 水 域 建设 用 地 未 利用 地 
耕地 1636095.10 1545.85 27207.35 760.73 444.76 4255.33 1372554.45 13365.60 173880.19 4162.97 1243.06 21058.31 
林地 2632.72 177889.46 16516.46 63.00 17.01 3761.44 40523.56 172234.56 22650.12 243.18 103.05 1551.31 
草地 30760.72 4130.78 1721320.01 1165.44 86.13 131127.51 223972.67 8299.49 1621368.70 3675.99 688.13 27335.35 
水 域 635.37 13.68 633.12 52136.74 26.64 55.53 3825.03 431.54 3054.37 42288.25 611.27 1080.97 
建设 用 地 ”1821.50 47.34 443.05 25.11 13724.81 295.99 24868.49 3763.56 45744.38 2120.55 13695.17 19428.62 
未 利用 地 4315.90 377.71 9664.04 575.97 13.50 449706.37 4621.07 2781.95 21815.02 366.65 17.73 394139.78 


候 变 化 和 土地 利用 之 间 存 在 复杂 关系 。 农 牧 交 错 
带 土 地 利用 类 型 对 气候 变化 的 响应 程度 不 同 , 且 有 
较 大 的 空间 异 质 性 ,林地 、 草 地 及 未 利用 地 受 年 平均 
气温 和 年 平均 降水 量 影响 较 大 ””。 自 然 环 境 因 子 变 
化 时 间 尺 度 较 长 ,对 短 时 期 内 的 土地 利用 类 型 变化 影 
响 较 小 约 。 榆 林 市 气候 呈现 暖 干 化 ,气温 降水 变化 不 
明显 9, 气候 变化 与 土地 利用 变化 关系 不 明显 。 
3.4.2 植被 因素 ”植被 覆盖 度 等 级 数据 与 土地 利用 
分 类 数据 属性 提取 至 样本 点 ,得 到 榆林 市 NDVI 与 
地 类 分 布 之 间 的 相关 关系 (图 4) ,耕地 与 中 高 覆盖 
度 的 相关 性 高 ,林地 与 中 覆盖 度 和 中 高 覆盖 度 相关 
性 较 高 ,草地 普遍 位 于 中 和 覆盖 度 及 以 下 等 级 ,高 覆 
盖 度 草地 ,未 利用 地 和 建设 用 地 与 低 履 盖 度 相关 性 
最 高 ,水 域 与 中 低 覆 盖 度 相关 性 高 。 植 被 覆盖 与 地 
类 变化 两 者 相互 影响 ,关系 复杂 。NDVI 数 据 经 计 
算得 到 植被 覆盖 度 等 级 (图 5$) ,地 表 禾 被 空间 上 呈 
现 均 衡 化 ,2000 年 后 , 低 .中 低 覆 盖 度 面积 减少 ,中 
覆盖 度 .中 高 覆盖 度 面 积 增加 ,高 覆盖 度 面 积 减 少 ， 
植被 覆盖 度 等 级 提高 。 


(a) 2000 年 À (b) 20104 


0 70 km 


RE ”中 低 PE PA 高 覆 
盖 度 AE EE TEE 
植被 覆盖 度 等 级 


图 4 2018 年 植被 覆盖 度 与 土地 利用 相关 关系 
Fig. 4 Correlation between vegetation cover 
and land use in 2018 


3.43 政策 因素 退耕 还 林 还 草 政策 直接 影响 土地 
利用 变化 ;部 分 容易 造成 水 土 流失 或 者 荒漠 化 的 耕 
地 停止 耕种 ,旱地 回 林 地 .草地 、 疏 林 回 林地 和 沙 地 
向 草地 的 转移 ? ;耕地 面积 以 2000 年 为 变化 节点 ， 
政策 实施 前 后 耕地 面积 先 增 长 后 急剧 减少 ,转移 
为 林地 和 草地 ;政策 实施 后 常用 耕地 面积 大 幅度 


À 


À (c) 2018 年 


70 km 


0 70 km 0 


ME (Rem CORE M 中 覆盖 度 国 中 高 覆盖 度 国 高 覆盖 度 


图 5 2000 一 2018 年 植被 覆盖 度 空间 分 布 图 
Fig. 5 Spatial distribution of vegetation cover during 2000—2018 
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减少 ,耕地 数量 和 质量 显现 不 稳定 性 ; 随 着 生 
活水 平 提高 ,对 供给 充足 .高 品质 和 无 污染 的 粮食 
需求 加 强 ,注重 协调 政策 驱动 下 耕地 数量 减少 与 耕 
地 质量 .生态 保护 之 间 的 关系 ”。 沙 地 治理 政策 与 
用 地 性 质变 化 直接 相关 ; 近 40 a 未 利用 地 空间 格局 
缩减 ,变化 率 较 大 , 沙 地 治理 过 程 中 ,面积 减少 。 为 
防止 进一步 沙化 ,进行 人 工种 植 林 草 ,沙漠 被 林 草 
覆盖 ;人 为 改变 沙 地 性 质 ,成 效 显著 , 区 域 生态 环境 
改善 ,避免 了 沙 地 恶化 的 延伸 。 

3.4.4 人 口 、 经 济 因 素 人 口 快速 增长 引起 土地 利 
用 的 范围 和 强度 变化 ;城镇 建设 用 地 、 农 村 居民 点 、 
交通 道路 以 及 其 他 各 类 公共 设施 建设 一 定 程度 上 
发 生变 化 ;1980 一 2018 年 ,西北 部 各 县 区 人 口 增长 ， 
榆 阳 市 .神木 市 . 绥 德 县 变化 最 显著 (图 6); 结 合用 
地 空间 格局 和 变化 ,建设 用 地 与 人 口 增长 分 布 区 相 
关联 。 西 南部 耕地 面积 增加 ,与 该 区 域 人 口 的 增长 
相关 。 社 会 经 济 发 展 的 影响 体现 在 土地 利用 类 型 
多 样 化 .各 用 地 类 型 面积 趋 于 平衡 和 建设 用 地 急剧 
增加 55 ;城市 扩张 并 占用 大 量 农业 用 地 ,导致 耕地 
面积 快速 缩减 后 ;早期 农业 产业 耕地 需求 增加 , 林 
草地 向 耕地 转移 ; 随 着 工业 发 展 和 城镇 化 加 速 , 建 
设 用 地 大 幅 增加 ;第 一 产业 变化 影响 耕地 需求 量 ， 
榆 阳 区 增长 率 最 大 ;第 二 产业 比重 较 大 ,成 为 推动 
全 市 经 济 发 展 的 主导 力量 六; 工业 和 三 次 产业 产量 


ChinaXiv 合 作 期 刊 


] 空 间 格局 变化 及 影响 因素 分 析 


增长 指向 神木 市 . 榆 阳 区 和 府 谷 县 ,与 建设 用 地 变 
化 相关 。 


4 讨论 


榆林 市 位 于 生态 环境 较为 脆弱 的 农 牧 交错 带 ， 
土地 利用 空间 格局 关乎 土地 利用 结构 优化 配置 和 
生态 环境 的 有 效 治理 ,因此 ,探讨 生态 环境 问题 背 
后 伴随 土地 利用 的 变化 和 区 域内 部 差异 成 为 需 
要 。 首 先 ,榆林 市 生态 环境 质量 对 土地 利用 变化 较 
敏感 。 有 研究 表明 ,1985 一 2000 年 榆林 市 土地 利用 
变化 使 生态 价值 总 量 增加 ,土地 生态 风险 降低 "™， 
北部 风沙 滩地 区 土地 利用 变化 过 程 相 对 活跃 ”; 本 
文 研 究 发 现 1980 一 2000 年 林地 面积 小 幅度 增加 , 草 
地 变化 强度 最 大 ,北部 草地 密度 核心 西南 扩展 ,与 
前 人 研究 结果 相 一 致 。 其 次 ,土地 利用 分 布 差 异 导 
致 生态 效应 的 差异 性 。 榆 林 市 各 土地 利用 类 型 分 
布 格局 差异 明显 ,各 县 区 生态 风险 差异 显著 2” , 土 
地 综合 承载 力 区域 差 异 较 大 '。 最 后 ,榆林 市 人 为 
活动 成 为 土地 利用 变化 不 可 忽视 的 因素 ;前 人 研究 
表明 经 济 发 展 对 土地 利用 变化 的 影响 最 为 显著 2 。 
土地 利用 空间 格局 和 长 时 间 尺 度 的 演变 趋势 问题 
得 到 回答 ,但 是 对 于 土地 利用 变化 产生 的 后 果 没 有 
具体 分 析 和 展示 ,因此 ,以 土地 利用 变化 为 基础 的 
生态 风险 时 空格 局 演变 及 评估 成 为 下 一 步 研 究 的 


人 口 /104 人 
= 19804 
= 20004 
= 20184 


第 二 产业 /108 元 
= 20004 
= 2010 年 
= 2017 年 


@) 第 一 产业 府 谷 县 


第 一 产业 /10 元 
= 200048 
= 20104F 
= 20174F 
子 洲 县 
(d) 第 三 产业 府 谷 县 
神木 市 佳 县 
me FP W/05 
横山 区 Rag "2000F 
= 20104 
定 边 县 É = 201744 
靖 边 县 清润 县 


子 洲 县 


图 6 榆林 市 经 济 产 量变 化 


Fig. 6 Changes in economic output of Yulin City 
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结论 

本 研究 结合 榆林 市 地 域 特殊 性 ,利用 核 密度 分 
析 法 ,分析 了 榆林 市 近 40 a 的 土地 利用 空间 格局 变 
化 ,并 揭示 了 土地 利用 变化 的 主要 影响 因素 ,得 出 
以 下 主要 结论 : 

(1) 榆林 市 近 40 a 土地 利用 变化 以 耕地 面积 减 
少 、 建 设 用 地 大 幅度 增加 和 未 利用 地 减少 为 主 。 空 
间 格 局 演变 特征 为 :耕地 2 个 核心 高 密度 位 于 西南 
部 和 东南 部 ;东南 部 耕地 密度 呈 缩 减 趋势 ,西南 部 
呈现 相反 的 变化 趋势 。 林 地 密度 核心 位 于 东南 部 ， 
星 “ 多 核心 增长 "。 草 地 密度 核心 位 于 北部 和 西南 
部 ,以 北部 核心 为 主 ; 北 部 核心 向 南 扩展 ,等 级 增 
高 ,西南 部 地 区 密度 等 级 逆向 变化 。 建 设 用 地 密度 
核心 位 于 西北 部 6 个 县 区 ,变化 幅度 大 。 未 利用 地 
以 西北 部 毛乌素 沙漠 为 核心 。 

(2) 气温 和 降水 量 的 增加 ;植被 覆盖 等 级 和 面 
积 变化 显著 ;人 为 因素 是 榆林 市 土地 利用 变化 的 主 
导 因 素 。 耕 地 变化 与 退耕 还 林 还 草 政策 的 实施 显 
著 相 关 ; 建 设 用 地 变化 与 经 济 产量 增加 和 人 口 的 增 
长 相关 。 毛 乌 素 沙 地 治理 下 未 利用 地 大 幅度 减 
少 。 退 耕 还 林 还 草 减 缓 了 水 土 流失 ,改善 了 生态 环 
境 。 经 济 发 展 , 人 为 活动 愈加 剧烈 ,建设 用 地 大 幅 
度 增加 。 

(3) 以 生态 修复 和 保护 为 主 ,实现 科学 .绿色 、 
可 持续 的 高 质量 发 展 。 坚 持 山 水 林 田 湖 草 沙 综合 
治理 ,加 强 水 土 保持 生态 建设 。 未 来 应 该 防止 耕地 
浪费 ,协调 好 政策 实施 与 耕地 质量 和 生态 保护 之 间 
的 关系 ,优化 传统 农 牧 业 。 根 据 沙 地 环境 合理 有 效 
地 改善 沙 区 环境 ,局 部 以 自然 恢复 为 主 , 减 少 人 为 
干扰 。 
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Spatial pattern changes and influencing factors of land use 
in Yulin City in the past 40 years 


LI Yali™, YANG Fenli, YANG Lian’an'**, SHANG Xiaoqing'*“, 
HU Guigui™, JIA Lujing 


(1. Shaanxi Key Laboratory of Earth Surface System and Environmental Carrying Capacity, Northwest University, Xi’ an 710127, 
Shaanxi, China; 2. Institute of the Yellow River Research, Xi an 710127, Shaanxi, China; 3. College of Urban and Environmental 
Sciences, Northwest University, Xi’an 710127, Shaanxi, China; 4. Academy of Agriculture Sciences of Xianyang, Xianyang 
712000, Shaanxi, China) 


Abstract: Yulin City is located in the northern agro-pastoral zone of China and has a fragile ecological environment. 
It is an important barrier to ecological security. Land use patterns are affected by the allocation structure and 
management of the ecological environment. In this study, the land use patterns of Yulin from 1980 to 2018 were 
analyzed to identify the relationship between spatiotemporal changes and influencing factors. The primary data 
source was land use data at a 30-m resolution. The normalized difference vegetation index was used to represent 
vegetation. Data were taken from three industries to represent economic factors. Demographic data were used to 
reflect demographic changes. Kernel density analysis was used to calculate the density of data points with land 
use attributes. Different indices were used to analyze the speed, intensity, and structure of changes in land use 
patterns. The results showed that changes in land use patterns were mainly due to a decrease in cultivated land, 
substantial increase in construction, and decrease in unused land. Cultivated land decreased in the southeastern 
area but increased in the southwest. Forestland showed “multicore growth”. The grassland was most dense in the 
north and increased over time. Construction was most dense in six counties in the northwest and showed a large 
increase over time. Most unused land was in the northwest, which included the Mu Us Desert, and was 
significantly affected by land management. The spatial changes in land use patterns were attributed to the 
combined effects of policy, socioeconomic, and demographic factors. The results indicate that the waste of 
cultivated land should be prevented. Additionally, the relationship between policy implementation, cultivated land 
quality, and ecological protection should be coordinated, and traditional practices of agriculture and animal 
husbandry should be optimized. Man-made construction and development in desert areas should be reduced. The 
interweaving effects of human activities and changes in the ecological environment should be considered to 
promote ecological protection. We advocate the concept of ecological civilization, which can be defined as the 
mutual symbiosis between man, nature, and the geographical environment. 
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